Dislocations and the critical endpoint of the melting line of vortex line lattices 
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We develop a theory for dislocation-mediated structural transitions in the vortex lattice which 
allows for a unified description of phase transitions between the three phases, the elastic vortex 
glass, the amorphous vortex glass, and the vortex liquid, in terms of a free energy functional for 
the dislocation density. The origin of a critical endpoint of the melting line at high magnetic fields, 
which has been recently observed experimentally, is explained. 
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Since the pioneering work Q where the first-order flux 
line lattice (FLL) melting into an entangled vortex liquid 
(VL) was established, there has been a continuous devel- 
opment of our views of the vortex lattice phase diagram 
in high-Tc superconductors Weak point disorder was 
shown to drive the vortex lattice into a vortex glass (VG) 
state with zero linear resistivity Observations of 

hysteretic resistivity switching and magnetization mea- 
surements H have supported the first order melting of 
very clean lattices. A crossover from the first order melt- 
ing at low magnetic fields to a continuous VG-VL transi- 
tion has been related to the structural transition within 
the vortex solid which is described f^flO[ as a topologi- 
cal transition between the low-field elastic VG, free from 
topological defects pl] , p2| and maintaining quasi long- 
range translational order and the high field amor- 
phous VG, where disorder generates proliferation of dis- 
locations Jll| . A simple picture of the vortex phase dia- 
gram has emerged where the three generic phases - VL, 
the high field amorphous VG, and the low field, low tem- 
perature quasi-lattice or Bragg glass (BrG) |]l4|| - are gov- 
erned by the three basic energies: the energy of thermal 
fluctuations, pinning, and elastic energies. The transition 
lines are determined by matching of any of the two basic 
energies, and the match of all three energies marks the 
tricritical point where the first-order melting terminates 
P|. While this simplistic picture is supported by obser- 
vations on BSCCO, it fails to describe the YBCO phase 
diagram where the endpoint of the first-order melting line 
appears to be separated from the point where topological 
transition and melting line merge [15 . Furthermore, the 
Lindemann-criterion approach of [ ^~ W| can only quali- 
tatively estimate the location of the transitions telling 
nothing about the nature of the vortex states involved. 

In this Letter we undertake a step towards a quanti- 
tative theory of the entire vortex lattice phase diagram 
based on a unified description. We propose a model where 
all the transitions are dislocation-mediated, and derive the 
free energy for an ensemble of directed dislocations as 
a function of the dislocation density in the presence of 
disorder and thermal fluctuations. Each of the experi- 
mentally observed phases is characterized by its inherent 
dislocation density or, equivalently, by the characteristic 
dislocation spacing Rrj. The elastic VG is dislocation- free 



and has Rd = oo. The VL can be viewed as a vortex ar- 
ray saturated with dislocations such that Rd ~ a, and in 
the amorphous VG, Rd ~ R a , where R a is the so-called 
positional correlation length on which typical vortex dis- 
placements are of the order of the lattice spacing a [0. 
Each phase corresponds to one of the local minima in the 
dislocation ensemble free energy, and dislocation densi- 
ties in these minima represent the equilibrium dislocation 
densities in the corresponding phases. The global mini- 
mum corresponds to the thermodynamically stable phase 
under the given conditions, phase transitions occur when 
two local minima exchange their role as global minimum. 
This mechanism for the transitions enables us to derive 
Lindemann-criteria both for the locations of the thermal 
melting line and the disorder-induced instability line of 
the BrG. The characteristic scale set by the mean distance 
between free dislocations offers a natural explanation of 
the critical endpoint of the first-order melting line. While 
at low magnetic fields R a 3> a and the amorphous VG 
appears to contain significantly less dislocations than the 
VL, at higher field where R a — a the two phases become 
thermodynamically equivalent and the first order melting 
line terminates. 

A free energy for the dislocation degrees of freedom gov- 
erning phase transitions is derived from the vortex lattice 
elasticity theory. Dislocations in the FLL can be both of 
screw or edge type, but in cither case they are confined 
to the gliding plane spanned by their Burger's vector b 
and the magnetic field |L6| . The single dislocation energy 
consists of the core energy E c and of the logarithmically 
diverging contribution from the long-range elastic strains 
|l7| ] . Accordingly, the dislocation ensemble is modeled as 
an array of elastic lines with a long-range Coulomb-like 
interaction. The energy penalty for Burger's vectors with 
b > a leads to an effective hard-core repulsion of dislo- 
cations which conspires with the planarity constraint to 
limit the entropy gain of dislocations. This favors directed 
dislocations threading the entire sample [Q rather than 
large dislocation loops |is[| . A single directed dislocation 
line is parameterized by its displacement field ud(z) and 
described by the Hamiltonian 



H D [u D ] = j dz (e s + ^e D (d z u D ) 2 ^ 
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where the stiffness to oc In (l/k z a) has a logarithmic dis- 
persion from the long-range interaction and E s is the 
self-energy of a straight dislocation. For thermal melt- 
ing, the short wavelength limit, k z ~ 1/a, is relevant. 
After rescaling z — z\\/ c^/cqq dislocation energies be- 
come isotropic with er> = E c = cpKb 2 /4ir (cd rs 1) and 
E s = E c + Kb 2 /4:TT\n(L/a) where K = ^c^cee is the 
isotropized elastic constant (044 and C66 are the tilt and 
shear moduli of the vortex lattice, respectively. L is the 
system size in the transverse direction). Note that the 
Peierls barrier W p and the associated "kinking" |l7j of 
dislocation lines can be neglected near the melting tran- 
sition. It can be shown that kinks are irrelevant above the 
temperature Tj. ~ a^/Wd^d- Since we find W p < 10~ 4 i? c , 
Tfc is much lower than T m . Therefore, the basic length 
scale along the magnetic field is solely set by the compe- 
tition of FL tilt and shear and given by a z ~ yj C44/ 1 cqq/2 
(a z w a in the rescaled system). The free energy of a 
single dislocation can be readily calculated from the par- 
tition sum Zd = J Dud exp (— (3Hd[ud]) by Gaussian 
functional integration and consists of the core energy, the 
long-range strains elastic energy, and the entropy term: 

F D (L) Kb 2 , f L\ 1 , / 2nTa z \ , , 
- E c + — In ( - ] - T- — In [ 1 + ± ) (2) 
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A second order phase transition due to the formation of a 
single dislocation is prohibited by the logarithmic diver- 
gence of its elastic energy. However, in an ensemble of dis- 
locations, this divergence is screened out on the distance 
l/2pa, where 2p is the dislocation density (the ensem- 
ble is topologically neutral to avoid the accumulation of 
stresses). An additional entropy cost (oc p 3 ) comes from 
the steric repulsion, and the resulting free energy density 
then reads 



f(j>) = 2p (E c 
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f(p) can be derived in a more rigorous manner by map- 
ping dislocations onto a quantum system of 2D Fermions 
with Coulomb interaction pO| , pl[ . At high dislocation 
densities, the screening mechanism in (^) is purely en- 
tropic in nature and equivalent to the formation of an 
exchange-hole in the 2D fermionic system. A first order 
melting following from (0) occurs at T m « 1.5E c a z ~ 
0.15Ka 3 , which is equivalent to melting according to the 
Lindcmann-criterion with a Lindemann-number cl ~ 0.2; 
in good agreement with experimental and numerical re- 
sults. At the melting transition dislocations proliferate 
with a high density p m w 0.3a -2 , hence the VL is satu- 
rated with dislocations. The presented scenario does not 
require a simultaneous proliferation of disclinations, how- 
ever, the latter are likely to appear at the high dislocation 
densities involved. 

In the presence of a random potential V P i n (r), the 
collectively pinned dislocation-free vortex array passes 



through three different scaling regimes: (i) Small scales 
where vortex displacements u are smaller than the co- 
herence length £ and perturbation theory applies p3[ . 
(ii) Intermediate scales where £ < u < a and disor- 
der potentials seen by different FLs are effectively un- 
corrected. This regime is captured in so-called ran- 
dom manifold (RM) models pjl4|, leading to a roughness 
G(r) = ((u(r) - u(0)) 2 ) « 4(a/27r) 2 (r / ' R a ) 2C - RM where 
( RM w 1/5 for the d — 3 dimensional RM with two 
displacement components. The crossover scale to the 
asymptotic behaviour is the positional correlation length 
R a where the average displacement is of the order of the 
FL spacing: u ~ a/2ii |p2f . (iii) The asymptotic Bragg 
glass regime where the et-periodicity of the FL array be- 
comes important for the coupling to the disorder and the 
array is effectively subject to a periodic pinning poten- 
tial with period a flg)| . Here the logarithmic roughness 
G(r) « 4(a/27r) 2 ln(er/i? a ), i.e., ( BrG = 0(log) ||Jl4| 
takes over. 

In a disordered system at T = the mechanism for dis- 
location proliferation is fundamentally different from the 
thermal melting discussed before. While thermal melt- 
ing is driven by the entropy gain from unbinding dislo- 
cations, at T = the FLL benefits energetically from 
adjusting itself to the disorder, and dislocation prolifera- 
tion leads to disorder energy gain through the dislocation 
degrees of freedom. It has been shown in Refs. [|ll|jl2| 
that the BrG phase in 3D is stable with respect to dis- 
location formation. Instabilities arise, however, from the 
sub- asymptotic regimes [in Ref. [jTl| this has been partly 
taken into account by considering displacements in the 
dislocation core]. To handle analytical difficulties and to 
provide a unified treatment through all scaling regimes, 
we develop an approach to the 3D problem based on an ef- 
fective random stress model which has the same displace- 
ment correlations as the full non-linear disordered model 
but allows for a separation of dislocation and elastic de- 
grees of freedom. This idea is motivated by the renormal- 
ization group (RG) for the 2D BrG which explicitly shows 
it renormalizes asymptotically into a random stress model 
pa] and has been used in Ref. 24 to show the instability 



of the 2D BrG with respect to dislocations. For simplicity 
we consider a uniaxial model (in the incompressible limit 
Cn ^> egg) which yields the same dislocation energetics as 
the isotropized two-component model. The Hamiltonian 



H[u] 



-K(Vu) 2 + a ■ Vu 



(4) 



where <r(r) is the random stress field which we as- 
sume to be Gaussian distributed with a second moment 
(Ti(k)o-j(k') = (%£(fc)(27r) 3 <5(k + k') characterized by the 
function £(&;) in Fourier space. The effective random 
stresses causing displacements with the same (2-point) 
correlations as those for the RM or BrG regime are 
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BrG: ^k^a 2 

RM: B 1/5 K 2 k~ 1 a 2 (kR a ) 
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with a numerical constant B 1 / 5 depending only on the 
roughness (rm = 1/5 (the exact crossover between the 
two regimes is non-trivial p4| ). The validity of the ran- 
dom stress model is well-established in 2D. Besides, the 
functional RG treatment of the BrG in d — 4 — e di- 
mensions shows that displacements asymptotically obey 
Gaussian statistics up to the first order in e [ p7| , which 
can always be modeled by an effective random stress field. 
Similarly, in a real-space RG analysis pq |, the relevant 
fixed point is perturbative in e and an analogous argu- 
ment applies. 

Starting from (^) the energy for an ensemble of curved 
dislocation lines with a Burger's vector density b(r) in the 
disordered system can be calculated (for example along 
the lines of Ref. ^5|). In the random stress model the 
Hamiltonian decouples into the elastic part and a disloca- 
tion part: 



H D [b] =J J t y b r ' br'G (r - r') + Jh r ■ g r 



(6) 



where Gq(t) — l/(47rr) is the 3D Green's function. Here 
g(r) is an effective random potential for dislocations lines 
defined by the transversal part of a through V x g = gt 
(cf. Q]). This energy contains the long-range elastic 
energy E s of dislocations in the first term and in the 
stochastic second term the disorder energy E^is of the 
dislocation which allows dislocations to gain energy by 
optimizing their paths. The disorder energy is completely 
determined by the FL displacement correlations through 
.g^k^k') = <5y J ft:S(fc)fc- 2 (27r) 3 (5(k + k / ) in the different 
regimes given by (ph. For a single directed dislocation line, 
the Hamiltonian (Q) reduces to the problem of a directed 
elastic line (with a logarithmic dispersion, see above) in 
a random potential that is long-range correlated due to 
@ . For a directed dislocation of length L z and confined 
in the transversal direction to a scale L, the elastic en- 
ergy is again given by (g) with T = 0. The mean square 
disorder energy fluctuations are 



E 2 dls (L Zl L) = 



BrG: cbtgE 2 d .L z L 
RM: c RM E 2 D L z L^ 
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(7) 



with numbers cbtG, c rm = These expressions give 

an estimate of the typical disorder energy a dislocation 
line can gain. They neglect rare fluctuations in the en- 
ergy gain from optimally positioning the dislocation in 
the transversal plane which give logarithmic corrections ~ 
O(lnL) |p^| . The optimal path of the dislocation will be 
rough up ~ L^ Z D with an exponent C,d ■ The roughness can 
be obtained by a Flory argument that equates the elastic 
energy from the deformation ££>(L)u 2 D /L (with a loga- 
rithmically dispersive ejj ~ lnL on large scales) and the 

cV2 



scale set by the fluctuation wavelength L. This yields 
( D (BrG) = 1 - Oilog 1 ' 2 ) and ( D (RM) = ^ ~ 0{log 1 / 2 ) 
where logarithmic corrections come from the dispersion of 
the stiffness and rare fluctuations. Since (fj(BrG) < 1, 
the BrG appears to be marginally stable against penetra- 
tion of a single directed dislocation whereas (,d{RM) > 1 
such that the random manifold is clearly unstable. Note 
that the scaling arguments of Ref. taking into ac- 
count rare fluctuations give the same result as the sim- 
ple Flory argument neglecting rare fluctuations. We ar- 
gue that the sub-asymptotic instability of the FLL in 
the RM regime on scales L < R a , combined with the 
asymptotic stability in the BrG regime for scales L > R ai 
leads to a disorder-induced dislocation proliferation via 
a weak first order phase transition. The characteristic 
dislocation density p c ~ R~ 2 at the transition is given 
just by the crossover scale R a - The discontinuities in 
this transition are small and may eventually disappear 
if the length scale R a becomes of the order of typical 
sample dimensions. Note also that, qualitatively, this re- 
sult is based only on the fact that the instability sets 
in within a sub-asymptotic regime; the random stress 
model is only used to quantify our findings. It enables 
us to estimate typical free energy minima for ensembles 
of dislocation lines with rough optimized paths at T = 0. 
The screened long-range elastic energy density for a (neu- 
tral) dislocation ensemble with density 2p is given by 
e D (p) = 2p(E D + (Kb 2 /4ir) In {l/y/2ap 1 ' 2 )) as in (g) at 
T = 0. Dislocations are confined to a transversal scale 



Rd 
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set by the distance to the next dislocation, 



hence they optimize their disorder energy gain on each 
longitudinal scale L z ~ Ru independently. Using (j^) with 
L z = L = p~ x l 2 for the BrG regime (p < R^ 2 ) and the 
RM regime (p > /2~ 2 ), we can estimate the corresponding 
minimal free energy densities (not considering logarithmic 
corrections from rare fluctuations) 



f(p) ~ e D (p) 



BrG: 2A BrG E c p 
RM: 2A RM § (pa 2 )% 



(8) 



with numbers AbtG^Arm = When both results 

in (||) are combined one indeed finds a local minimum 
in the free energy density at p ~ R~ 2 that characterizes 
an amorphous VG phase. Over a wide range of mag- 
netic fields the dislocation density in the amorphous VG 
is much lower than in the VL for which we have found 
p w 0.3a~ 2 above. The elastic BrG phase looses stabil- 
ity with respect to dislocation proliferation and a tran- 
sition into an amorphous VG phase if the local mini- 
mum at p R~ 2 becomes the global free energy mini- 
mum. This occurs via a weak first order transition when 
Ra/o, — C with a "Lindemann-number" C = O(10). 
This is identical to the Lindemann criterium obtained in 
Rcfs. 11 HJ] including the numerical value of C j22|. As 
it was shown in Ref. pOl it is equivalent to a Lindemann- 



criterion ((u(a) — it(0)) 2 



disorder energy E 2 is (L z = L,L) on one large length (see pfl) with cl — j^C 1 / 5 w 0.1. 



c\a 2 in its more familiar form 
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FIG. 1. Schematic phase diagram of YBCO. Insets show 
typical free energy densities / of a dislocation ensemble 
as function of the dislocation density p. 

We have identified three possible minima for the free 
energies (||) and (||): (i) The dislocation-free minimum 
at p = which is stable in the elastic BrG phase at low 
T and low H. (ii) The minimum at p ~ a~ 2 that be- 
comes stable in the disorder-free case for high T in the 
VL. (iii) A minimum at p ~ R~ 2 which is realized in the 
amorphous VG. Combining our results for the thermal 
melting and the disorder-induced "melting" , we have ob- 
tained a qualitative theory for the entire phase diagram 
of the vortex matter. The experimentally observed criti- 
cal endpoint of the first order melting line finds a natural 
explanation. At elevated fields the positional correlation 
length R a decreases fl(i[| and finally reaches R a ~ a such 
that the two free energy minima of the VL and the amor- 
phous VG merge. Both these phases become thermody- 
namically indistinguishable and have identical equilibrium 
lattice order. Above the critical endpoint there might 
still exist a dynamic transition (or crossover) which in- 
volves the thermal depinning of dislocations, similar to 
the well-known thermal depinning transition of, for ex- 
ample, a single pinned vortex line. The critical endpoint 
is located at the magnetic induction determined by the 
condition R a « 2a (which is again equivalent to a Linde- 
mann criterion but with C cp ~ 2 or cl ~ 0.14). This gives 
the estimate B cp /H c2 » (2tt)- 13 / 3 (5/e)- W/d C~ 16/15 ~ 
O(10 _1 ) where we followed Ref. || and introduced a di- 
mensionless disorder strength S/e ~ 10~ 3 for high-T c ma- 
terials (cf. jL0|). This value is approximately a factor 10 
higher than the instability field BbtG of the BrG and thus 
the "coexistence point" where the topological transition 
line ends in the first order melting line and all three phases 
- elastic BrG, amorphous VG, and VL - can coexist. 
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